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Abstract. We study the nonmesonic weak decay (NMWD) AN — > nN of the ]^C 
hypernucleus induced by the nucleon N = n,p with transition rate Fjv. The nuclear process 
fj ' is described by the interplay of two models; one describing the NMWD of hyperon A in the 

r^ , nuclear environment, and the other taking into account the Final State Interaction (FSI) of 

r^ ' the two outgoing nucleons with the residual nucleus. The first one is done in the framework 

of the Independent-Particle Shell- Model (IPSM), with the decay dynamics represented by the 

exchange of tt+tj+K+p+oj+K* mesons with usual parametrization. For the second one is used 

a time dependent multicoUisional intranuclear cascade schema (implemented in the CRISP code 

^1. , - Collaboration Rio-Sao Paulo). The results obtained for inclusive and exclusive kinetic energy 

I ■ spectra, and the angular correlation are compared with recent data from KEK and FINUDA 

fvq ' experiments. The calculated ratio (Fn/Fp)^ , between the numbers of emitted back-to-back 

fT^ , nn and np pairs, is in good agreement with the experimental data. 
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1. Introduction 

Since tlie mesonic channel A -^ ttN is strongly inhibited inside nuclei due to Pauli blocking, the 
K> ■ nonmesonic weak decay (NMWD) AA^ -^ nN emerges as the dominant channel. This primary 

j^ I decay can be induced by a neutron (An -^ nn) or a proton {Ap — )• np) with widths r„ and Tp, 

respectively. During the last two decades quite significant efforts have been invested in solving 
the so called 'NMWD puzzle'. This puzzle is related to the large discrepancy between earlier 
experimental data and theoretical predictions on the NMWD spectra, that yield information 
about both the primary nonmesonic decay, and the subsequent Final State Interactions (FSI) 
originated by the nuclear medium. The main aim of the present work is to analyze in which 
way the FSI modify the spectra obtained within the Independent Particle Shell Model (IPSM), 
and to find out to what extent it is possible in this way to explain the experimental data. The 
spectra that we are interested in are: 1) the single-nucleon spectra S]\f{EN), as a function of 
one-nucleon energy S^r, and 2) the two-particle-coincidence spectra as a function of: i) the sum 
of kinetic energies En + E]y = EnN, S^iEnN), and ii) the opening angle 9nN, Sj\ficos9nN)- The 
explicit relationships for these transition probability densities can be obtained by performing 



derivatives on the appropriate equation for Tj^ [H El El H] , 
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The measured spectra are obtained by counting the number of emitted nucleons AN^r within the 
energy bin AE = 10 MeV, or the angular bin Acos0 = 0.05, always corrected by the detection 
efficiency. Here we take advantage of the fact that in most of the KEK experiments O El [7] 
AN^r are normahzed to the number of NMWD processes, NtvM; while for the FINUDA proton 
spectra [3 |9] we have at our disposal also the number of produced hypernuclei Nj^. Therefore, 
we can explore the following relationships 
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where AFjv is the emission rate of protons within the experimentally fixed bin, and Ti\im and 
Th are, respectively, the NMWD rate and the total hypernuclear decay rate. It could be worth 
noting that Th = Tm + Tnm, i-e., is the sum of mesonic and nonmesonic decay rates. It is not 
clear whether one should use experimental or theoretical values for the rates Tjvm, and Th, and 
we will adopt the first possibility. Then the correspondences between the theory and data are 
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Tnm = 0.95 ± 0.04 [6], and Th = 1.21 ± 0.21 [IQj. The theoretical decay widths 
^ are: ATn{En) = SNiEN)AE, ATiyiEnN) = 5'Ar(£'„Ar)A£', and ArAr(cos6'„Ar) = 
S'Ar(cos^„Ar)Acos0 for the spectra explained in 1) and 2). As the one-proton (one-neutron) 
induced decay prompts the emission of an np (nn) pair the total neutron kinetic energy width is 
Ar„t(^„) = {Sp{En) + 2Sn{En))AE. For ATp{Ep), ATnt{En), and Ar^(cos0„iv) are available 
the data, with the corresponding errors. As for the kinetic energy sums, only the data for 
yields AY^iEnN) were reported so far, without specifying their errors [5J. It is not known how 
ATj\f{EnN) and AYj\f{EnN) are related with each other, and thus one is forced here to normalize 
the theoretical results to the data in a similar way as done in Ref. [2], i.e., 
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where 



Y^-P = Y^AY^^P{El,^), f^ = Y^SN{El,^)AE, 



(4) 
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with m being the number of kinetic energy sum bins, and the barred quantities indicate that 
and they are constructed with events constrained to E^ > 30 MeV, and cos OnN < —0.7. 



2. Independent Particle Shell Model for the primary NMWD 

Within the IPSM the spectra defined in (P) read [H El El S]: 
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where the summation goes over aU occupied single-particle states j = {n/j} for a given N = p,n, 
and 



E'^ = {A - 2){A - l)Aj^ - En[{A - if - cos^ 9nN], 



(9) 
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A is the nuclear mass, e's are the single-particle energies, and A = Ma — M is the mass difference 
between the A-hyperon and the nucleon. The quantity Fj^^ {pP) contains the isospin dependence, 
as well as the dynamics, which is described by exchanges of complete pseudoscalar (tt, K, rj) and 
vector {p,uj,K*) meson octets with the weak coupling constants from Refs. |1HI12J. It is defined 
in Eq. (42) of Ref. ^, and depends on center of mass and relative momenta of the nN pairs. 
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with pm and pN being the momenta of the emitted particles. 



3. Monte Carlo calculation of FSI with CRISP 
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Figure 1. Two-dimensional spectrum for protons S(j)p,cos6np) as evaluated from the IPSM 
Eq. (USD. 



After having modelled the primary process, we investigate the FSI through the time evolution 
of Monte Carlo samples constructed to represent the decaying nuclear system, making use of a 
time dependent multicollisional intranuclear cascade schema (implemented in the CRISP code - 
Collaboration Rio-Sao Paulo). The initial sample configurations are spherically shaped portions 



of cold Fermi gas with A — 2 nucleons, where are incorporated the two energetic intruder nucleons 
resulting from the primary process, with the IPSM spectra 

Sn{pn, cosOnN) = -, — 1 2 — = (-4 - 2) — y^PNPl^JNipP)^ (13) 

dp]\jd COS OnN 71" 

which are shown in Fig. [1] for the emission of the np pair, and are the starting configuration 
for the time evolution of the 74-nucleon system, following the sequence of all possible binaries 
collisions between nucleons. A square well potential is included to mimic the nuclear surface, 
keeping the low energy nucleons bound within the system. The depth of this potential well 
is properly chosen in order to have an unperturbed initial configuration without nucleon loss. 
In fact, the potential confines low energy nucleons through internal reflections, and permits 
refractions/deflections of energetic nucleons, as dictated by the tunnelling probability of the 
charged particles through the nuclear Coulomb barrier |13| . 

The time order of the nucleon-nucleon collisions during the evolution of the system is 
reconstructed after each collision between pairs of particles, and their displacement inside the 
nucleus. Consequently, the many-body configurations are followed step by step in phase space, 
and in time. The cascade phase is ended when no particle is able to leave the nuclear volume, 
and the number of bound nucleons in the residual nucleus does not change anymore. However, 
the system still can present some excitation energy compared to a cold nucleon Fermi gas. A 
thermal decay chain is performed to cool down the post cascade nuclear system by emission 
of protons, neutrons, and alpha particles, in competition with the fission process |14j . Briefly, 
the FSI were calculated through a Monte Carlo simulation of all nuclear processes after the 
primordial NMWD [13^ll4 1 [T5 l ll6 t [T7]. including important features, such as the Pauli blocking 
within both elastic and inelastic nucleon-nucleon collisions [13lll7j. 

Note that the interplay between the IPSM and the intranuclear cascade calculation is done 
through the initial configuration for the Monte Carlo simulation, since the two intruder nucleons 
in initial sample configurations have their energies, momenta, and opening angle determined 
stochastically from Eq. ()13p . Another important bridge between the two models is the choice of 
the isospins for the two nucleons that start the cascade process. This is also done via a Monte 
Carlo sampling of the outgoing channel according to the IPSM value for the Tn/Tp ratio in the 
primary NMWD process. 

4. Result and discussion 

In Fig. [2] are compared the calculated inclusive proton (left panel) , and neutron (right panel) 
spectra with the KEK [6l[7], and FINUDA [8, 9j data. It is well known that the SFI modify the 
corresponding IPSM spectra Arj\f{E]\f) = S]\f{Ej\f)AE, with S^iEN) given by ([6]), by promoting 
particles from the high energy region towards the low energy region. From the figure one can 
see that the FSI mechanism is engendered by the interplay between contributions coming from 
the cascade regime emission, and those induced by the evaporation of the warm residual post- 
cascade nucleus. The enhancement at low energy induced by the first effect agrees with data, 
while the quenching at high energy produced by the second effect is too strong, in disagreement 
with the data. We also note that the threshold for the opening of the cascade effect is sharper 
in the pn that in the nn case due to the Coulomb barrier. 

In left and right panels of Fig. [3] we show the angular distributions for the nn and np 
channels, respectively. We get reasonable agreement with experimental results in the case of 
the nn channel, but our calculation does not reproduce the peak observed experimentally in the 
back-to-back region for the np channel. This is due to too strong FSI quenching on the IPSM 
spectra Sp{cos6np) given by ([7]). 
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Figure 2. Calculated inclusive single particle spectra (full histogram), compared to the data 
from Ref. [U", 7] in open circles, and from Ref. [SJ E] in full circles. The dotted histogram are 
the contribution coming from cascade regime of the process and the dashed one depicts the 
contribution of particles evaporation of the warm post-cascade residual nucleus. 
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Figure 3. Angular distribution of the np-pairs (left panel), and the nn-pair (right panel) 
compared to experimental data from Ref. [6j[7]. 



In the right panel of Fig. H] we show the total kinetic energy of the correlated np pairs, and 
compared them to the KEK experimental data [5]. The theoretical predictions are normalized 
to the data according to Eqs. (JH), and dS]). We observed that both the centroid and width 
of the theoretical spectrum are quite similar to the experimental one, covering a broad energy 
interval from ~ 80 MeV up to the Q- value at ~ 150 MeV. For the nn channel, as shown in the 
left panel of Fig. HI the agreement between calculated and data is not so good as in the the np 
case. However, we would like to draw attention to the fact that, because of very low statistic, 
which in turn is due to difficulties in detecting two neutrons in coincidence, the experimental 
histogram contains relatively small number of events. 



We have also calculated the ratio between r„ and Tp by considering the emitted nn and np 
pairs in the back-to-back opening angle region cos 6nN < —0.7, getting the value (Tn/TpY^^ = 
0.454, in nice agreement with the recently reported experimental result (r„/rp)'^^P = 0.53 it 
0.13 [5j, where the background uncertainty is 0.05. At this point it is interesting to mention 
that the IPSM model calculation yields the value {Tn/TpY^^^ = 0.262. We can attribute the 
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Figure 4. The kinetic energy distribution of the nN pairs calculated for events with strong 
back-to-back correlation of the opening angles. The full histograms are the results from the 
calculation, and dotted histograms are data from Ref. ^. To define the back-to-back events we 
have adopted the same low energy cutoff, and the same angle interval as done in the experiment. 

difference between (r^/Fp)^''^^, and (r„/rp)^^' to a possible depletion in the proton emission due 
the Coulomb barrier at the residual nucleus surface, implying that a neutron produced during 
the primary process have higher probability of leaving the nucleus than a primary proton. 

5. Conclusions 

The role played by the FSI on the NMWD has been investigated. For the inclusive nucleon 
spectra we have depicted separately the contributions coming from the cascade phase, and the 
low energy ones coming from the particle evaporation process in the post cascade phase regime. 
The importance of the latter nucleons in reproducing the experimental data was put in evidence. 
In the analysis of angular distributions, we have obtained a reasonable agreement for the nn 
channel, while the results for the np channel still need to be improved. From the comparison of 
the F„/Fp ratios before and after the inclusion of FSI, one sees that these interactions increase 
the n/p ratio by a factor of ~ 1.8, yielding in this way the agreement with the measurement. 
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